arXiv:1503.03329vl [astro-ph.SR] 11 Mar 2015 


Astronomy Astrophysics manuscript no. ProtonatedPANH'Revised'Edited © ESO 2015 

March 12, 2015 


The electronic spectra of protonated PANH 

molecules 

J. A. Noble, C. Dedonder, and C. Jouvet 

CNRS, Aix-Marseille Universite, UMR-7345, Physique des Interactions loniques et 
Moleculaires (PIIM): 13397 Marseille Cedex 20, France, christophe.jouvet@univ-amu.fr 

Received March 12, 2015/ Accepted March 12, 2015 

ABSTRACT 

Aims. This study was designed to examine the viability of protonated nitrogen- 
snbstituted polycyclic aromatic hydrocarbons (IT''PANHs) as candidates for the carri¬ 
ers of the diffuse interstellar bands (DlBs). 

Methods. We obtained the electronic spectra of two protonated PANH cations, proto¬ 
nated acridine and phenanthridine, using parent ion photo-fragment spectroscopy and 
generated theoretical electronic spectra using ab initio calculations. 

Results. We show that the spectra of the two species studied here do not correspond 
to known DlBs. However, based on the general properties derived from the spectra of 
these small protonated nitrogen-substituted PAHs, we propose that larger H’''PANH 
cations represent good candidates for DIB carriers due to the expected positions of 
their electronic transitions in the UV-visible and their narrow spectral bands. 

Key words, astrochemistry - Molecular data - Ultraviolet: ISM - ISM: molecules - 
Methods: laboratory: molecular 


1. Introduction 


Polycyclic aromatic hydrocarbons (PAHs) are believed to be widely distributed th r oughout 


the i nterstellar medium (ISM) in both the gas and solid phases (jTielens 


2008; 


Salama . 


20081) . The observation of two sets of unidentified spectral features, the diffuse interstellar 
bands (DlBs) and the unidentified infrared bands (UIRs), is routinely attributed to the 
presence of PAHs (in absorption and emission, respectively). 

The DlBs - weak unidentified absorption features in the near ultraviolet (UV), visi¬ 
ble, and near infrare d (IR) region of the e lectromagnetic spectrum - were first observed 


in the 1920s (Heger 


1922 


Merrilll . 1193411 . These bands are observed on lines of sight 


containing sufficiently high column densities, such as those traversing diffuse interstel¬ 
lar clouds. Additionally, the UI R emission features at 3.3, 6.2, 7.7, 8.6, and 11.3 p,m, 
first observed in the 1970s (e.g. iGillett et al.l . Il973ll . are observed in the IR spectra of 


almost all objects. Although no individual carrier molecule has been identified with cer¬ 
tainty, the group of polyaromat i c hyd r ocarbons (PAHs) is believed to contribute to these 
bands (e.g. 


Allamandola et al 


1989; 


Salama et al 


19961 and references therein). The 
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develo pment of the so-called “PAH hypothesis”, which holds PAHs responsible for the 
UIRs (ILeger fc Pueetl . 1198411 due to UV-pumped IR fluorescence, was very quickly fol- 
lowed by an hypothesis suggesting that ionised PAHs were likel y carriers of the DIBs 


I Leger &: D’Hendecourt 


1985 


van der Zwet fc Allamandola 


1985). 


Both experimental and theoretical studies support the identificatio n of PAHs, and more 


speci f ically, PAH cations, as the carriers o f the DIBs and UIRs ( e .g. Brechignac,&_^inp , 


1998; 


2011 


Salama et al 


Gredel et al 


1998; 


Mattioda et al. 


2005 


Pathak fc Sarre 


200 . 


201111 . Models have shown tha t UV pumping (e .g. 


Mattioda et al. 


Garkusha et al 


Schntte et al 


1993) 


( 20081) 1. can drive MIR 


or even NIR absorption (in areas of low UV flux 
emission from PAHs, effectively linking DIBs with UIRs. One of the remaining points 
of debate is the nature of the PAH cations: whether they are radical cations (PAH’*') or 
protonated species (H+PAH). Radical cations have been shown experimentally to be highly 


react ive and thus will react readily with hydrogen to form protonated speci es (ISnow et al 


I998h : PAHs are, however, stripped of H atoms by interaction with photons (jTielens 


2008). 


More recent IR studies have suggested that the sub-group of nitrogen-substituted PAHs, 
or PANHs, exhibit spectral f eatures similar to those of PAHs, and may also contribute to 


unidentified spectral bands (IHndgins et al 


2008: 


Alvaro Galue et al. 


2005 


Bernstein et al. 


2005 


Mattioda et al 


20101) . One study in particular illustrates the importance of 


PANHs in photon-dominated regions (PDRs), as the spectra of PANH cations were not 
merely usef ul, but required , to fi t the 6.2 and 11 /im emission features observed towards 


NGC 7023 I Boersma et al 


2013h . This is not an unexpected result, as the abundance of 


N relative to C in the ISM is ^ 0.25 (|Spitzer 


1978 ) and thus PA NHs are likely relatively 


Hudgins et al 


(120051) conclude that the sub- 


abundant interstellar species. Additionally, 
stitution of N in PAHs could account for the variations observed in the peak position of 
the 6.2 p,m UIR. These spectroscopic studies in the IR have strengthened the case for the 
presence of charged PANHs in the ISM and their contribution to the UIRs, but very few 
studies have be en made of the elect ronic spectra of gas phase PANH cations for comparison 


with the DIBs (Drvza et al 


2012 ). 


PANHs have a high proton affinity, so it is highly likely that protonated PANHs exist 
in the ISM, particularly in ionised environments. It has been illustrated that, for a small 
sample of PANH molecules, the IR spectra of protonated PANH cations (H+PANH) better 
reproduce the 6.2 and 8.6 iiui UIRs than the spectra of ionised PANH radical cations 
(PANH+) (|Alvaro Galue et al.l . 120101) . However, to date no studies have determined the 
electronic spectra of H+PANH for comparison with the DIBs. 


In this study, we provide the electronic spectra of gas phase protonated cations acri¬ 
dine ([Ci 3 H 9 N]H+, hereafter AcH+) and phenanthridine ([Ci 3 H 9 N]H+, hereafter PhH+), 
nitrogen-substituted versions of the PAHs anthracene and phenanthrene (see Appendix lATl. 
We show that these protonated molecules present vibrationally resolved electronic states 
in the visible or near UV spectral region where DIBs are observed. 
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2. Methods 


2.1. Experimental methods 

The electronic spectra of the protonated aromatic PANH were obtained via parent ion 
photo-fragment spectroscopy in a c ryogenically - cooled quadrupole ion trap (Paul Trap 
from Jordan TOP Products, Inc.) (lAlata et al.l . 1201311 . T he setup is simi l ar to the one 
developed in several groups based on the original design bv lWang fc Wand ([20081). 

The pro tonated ions are produ ced in an electrospray ionisation source built at Aarhus 
University ([Anderson et al.l . 1200411 . At the exit of the capillary, ions are trapped in an 
octopole trap for 90 ms. They are extracted by applying a negative pulse of ca. 50 V and 
are further accelerated to 190 V by a second pulsed voltage just after the exit electrode. This 
timed sequence of pulsed voltages produces ion packets with durations of between 500 ns 
and 1 ps. The ions are driven by electrostatic lenses towards the Paul trap biased at 190 V 
so that the ions enter the trap gently, avoiding fragmentation induced by collisions. A mass 
gate placed at the entrance of the trap allows selection of the parent ion. The Paul trap is 
mounted on the cold head of a cryostat (Coolpak Oerlikon) connected to a water-cooled 
He compressor. Helium gas is injected in the trap using a pulsed valve (General Valve) 
trigge red 1 ms before the ions enter the trap, as previously reported by 


Kamrath et al 


( 201(1) . The ions are trapped and thermalised at a temperature between 20 and 50 K 
through collisions with the cold He buffer gas. The ions are held in the trap for several 
tens of ms before the photodissociation laser is triggered. This delay is necessary both to 
ensure thermalisation of ions and the efficient pumping of the buffer gas from the trap in 
order to avoid collision induced dissociation of the ions during subsequent extraction. 

The photo-dissociation laser is an optical parametric oscillator (OPO) laser, which has 
a 10 Hz repetition rate, 10 ns pulse width, a resolution of ^ 8 cm“^ and a scanning step of 
0.02 nm (laser from EKSPLA uab). The laser is shaped to a 1 mm^ spot to fit the entrance 
hole of the trap and the laser power is around 20 mW in the UV spectral region. After 
laser excitation, the ions are stored in the trap for a delay that can be varied between 20 
and 90 ms before extraction into the 1.5 m time-of flight mass spectrometer. The full mass 
spectrum is recorded on a micro channel plates (MCP) detector with a digitising storage 
oscilloscope interfaced to a PC. The photofragmentation yield detected on each fragment 
is normalised to the parent ion signal and the laser power. 


2.2. Calculations 


Ab initio calculation s were performed with the TURBOMOLE programme package 


I Ahlrichs et al 


19891), making use of the resolut ion-of-the-ide ntity (RI) approximation for 


the evaluation of the electron-repulsion integrals (IHattid . 120031) . The equilibrium geometries 
of the protonated species in their electronic ground states (So) were determined at the MP2 
(Mpller-Plesset second order perturbation theory) level. Adiabatic excitation energies of 
the lowest electronic excited singlet states (Si) we re deterniined a t the RI-ADC(2)(second 
order Algebraic Diagrammatic Construction level (ISchirmeil 119821) ). Calculations were per¬ 
formed with the correlation-consistent polarised valence double-zeta (cc-pVDZ) basis set 
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Excitation wavelength (nm) 

440 435 430 425 420 415 



T-'-1-'-1-'-1- 

22500 23OQ0 23500 , 24000 

Excitation energy (cm ) 



Fig. 1. Photofragmentation spectrum of protonated acridine compared to the simulated 
spectrum: a) the experimental spectrum plotted as the fragment ion signal (arbitrary units) 
as a function of the excitation energy; b) the simulated spectrum plotted as the Franck- 
Condon factors as a function of the excitation energy, calculated using the calculated 
ground and excited state frequencies. The main vibrational progressions (indicated in blue 
in panel a) involve, essentially, the two symmetric modes vio at 383 cm“^ (ring breathing) 
and Vis at 587 cm“^. These vibrational modes are illustrated schematically above the 
experimental spectrum, with arrows representing the direction and relative magnitude of 
the displacement of each atom. The intensities of the out of plane modes vi and are 
overestimated in the simulation. 


(|Woon fc Dunnind . Il993h . The vibrations in the ground and excited states were calculated 
for protonated acridi ne and phenanth ridine and the electronic spectra simulated using the 


PGOPHER software (Western 


201411 for Franck Condon analysis. 


3. Results and discussion 


3.1. Protonated acridine 

The protonated acridine photofragmentation spectrum is presented in Fig. [1^. The first 
electronic transition shows clear vibrational progressions in the visible starting at 442.24 nm 
(22612 cm“^/ 2.80 eV) and extending to ^ 400 nm, this absorp tion region being i n agree¬ 


ment with the absorption recorded in water solution at low pH (jRvan et al.l . Il997f) . These 
progressions involve the 383 cm“^ vibrational mode starting from the 0-0 transition (r'o), 
and a symmetric mode at 587 cm“^ (^^is)- The vibronic bands are narrow, with a FWHM 
of 10 ± 1 cm“^ that represents the convolution of the laser band width (~ 8 cm“^) with a 
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low temperature rotational contour i.e. the apparent line broadening due to the fact that 
the spacing between rotational lines is narrow compared to the laser width. 

The main fragmentation channels correspo nd to loss of H, 2 H/H 2 and m/z 28 (proba¬ 


bly corresponding to loss of H+HCN/H+HNC, 


Johansson et al 


(|201lll l complemented by 


weaker fragmentation channels: loss of m/z 26, 27, 29 and loss of m/z 51, 52 and 53 (see 
mass spectra in Appendix IbI) . 

Ground and first excited state calculations were performed, both with DFT/TD-DFT 
(Density Functional Theory/ Time Dependent-DFT) and MP2/ADC(2) methods, for the 
isomer protonated on the nitrogen atom, which is by far the most stable. The ground state 
energies of the isomers protonated on carbon atoms were also calculated using DFT with 
the B3LYP functional and the cc-pVDZ basis set, and they all lie more than 1.5 eV higher 
than the isomer protonated on N. The molecule stays planar in both states. The first excited 
state is calculated vertically at 3.13 eV, in ag reement with prev i ous C ASPT2 (complete 
active space perturbation theory) calculations (|Rubio-Pons et al.l . l200lh . and corresponds 
to a TTTT* HOMO-LUMO (highest occupied molecular orbital-lowest unoccupied molecular 
orbital) transition (^Lq in Platt’s notation). Optimisation leads to an adiabatic transition 
at 2.75 eV (2.60 eV when the difference in zero point energy is included), in agreement with 
the experiment (2.80 eV), which corroborates the location of the proton on the nitrogen 
atom. 

The changes in geometry between the ground and excited states are minor: for the outer 
rings, there is an alternation of long/short C-C bonds in Sq, which becomes short/long in 
Si, with a maximum difference in bond length of 0.053 A for the C1C2 and C7C8 bonds 
(see Appendix lAl for the So and Si optimised geometries and the numbering of atoms); for 
the center ring bearing the N atom, the NH bond is unchanged, NlOCll & N10C14 increase 
by 0.031 A, C9C12 & C9C13 increase by 0.035 A, while C11C12 & C13C14 decrease by 
0.008 A. The angles differ by at most 2°. 

The spectrum was simulated (see Fig. using the ground and exc ited state frequ encies 
calculated and the PGOPHER software for Franck Condon analysis (jWesternl . l2014^ . The 
agreement between the experimental and simulated spectra is quite good for the band 
positions, and allows assignment of the main vibrational progressions to the symmetric 
breathing mode vio (contraction/elongation of the rings along the long axis) alone or in 
combination with another symmetric mode I'ls (elongation of the center ring along the 
short axis with contraction of the outer rings). The scheme of the vibrations is shown at 
the top of Fig.[T] Two out of plane modes are weakly active in the experimental spectrum 
and overestimated in the simulation: vi is the butterfly mode and 1^3 mostly an out of plane 
motion of the NIO-H and C9-H atoms (seesawing motion). 


3.2. Protonated phenanthridine 

The protonated phenanthridine photofragmentation spectrum (presented in the upper 
panel of Fig. HD shows clear vibrational progressions in the near UV spectral region, with 
a first transition at 397.28 nm (25171 cm“^/ 3.12 eV). This transition is blue shifted 
by 2500 cm“^ compared to the AcH+ transition. Vibrational progressions on a mode of 
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Excitation wavelength (nm) 

400 395 390 385 380 



Fig. 2. Photofragmentation spectrum of protonated phenanthridine compared to the simu¬ 
lated spectrum: a) the experimental spectrum plotted as the fragment ion signal (arbitrary 
units) as a function of the excitation energy; b) the simulated spectrum plotted as the 
Franck-Condon factors as a function of the excitation energy, calculated using the calcu¬ 
lated ground and excited state frequencies. The main vibrational progressions (indicated 
in red in panel a) 1^5 and 1 ^ 17 , are illustrated schematically below the simulated spectrum. 


240 cm“^ starting on the 0-0 transition (t'o), and on a symmetric mode at 533 cm“^ (^ 17)1 
are observed for several quanta, and the bands are narrow, with a FWHM of 12 ± 1 cm“^ 
(see Appendix ICl) . As for AcH+, the main fragmentation channels correspond to loss of H, 
2 H/H 2 and H+HCN/H+HNC (see mass spectrum in Appendix IbI). 

Ground and first excited state calculations performed for the isomer protonated on the 
nitrogen atom show that PhH+ stays planar in both states (see Appendix lAl). As for AcH+, 
this isomer is by far the most stable, the calculated ground state energies of the isomers 
protonated on carbon atoms (DFT/ B3LYP/ cc-pVDZ) all lie more than 1.6 eV higher 
than the isomer protonated on N. The first excited state is a tttt* state calculated vertically 
at 3.45 eV and optimisation leads to an adiabatic transition at 3.22 eV, in agreement with 
the experiment (3.12 ev). 


3.3. Three-ringed protonated PANHs 

For both protonated acridine and phenanthridine, the proton is located on the nitrogen 
atom, as shown by the excellent agreement between experiments and calculations: the 
calculated excited state transitions are within 0.2 eV from the experimental band origins 
and the spectra simulated with the calculated frequencies reproduce well the experimental 
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Table 1. Transition origin (eV) for selected neutral and protonated PAHs and PANHs. 


PAH 

Transition 

PANH 

Transition 


origin (eV) 


origin (eV) 

Anthracene 

3.43^ 

Acridine 

3.22^’® 

AnthraceneH*^ 

2.52^^ 

AcridineH*^ 

2.80® 

Phenanthrene 

3.63® 

Phenanthridine 

3.64’’ 

PhenanthreneH*^ 

2.08®’® 

PhenanthridineH^ 

3.12® 


Re ferences. ( 1) 


(4) 

( 8 ) 


Alata et al. 


digsdl : (2) 


Alatal (l2012h : (9) 


Lar^ert et al 
( 2010b[l: (5) this work: (6) 


Prochorow et al 


Garkusha et al. 


Hager fc Wallace 

Il20ll 


1998^; (3) 


Rubio-Pons et al 


( 1988l ~l: (7) llProchorow et al 


( 2001 ); 


200J); 


progressions. Protonated acridine and phenanthridine are planar in the ground and excited 


states, which was not the case for protonated anthracene and phenanthrene (|Alata et al 


2010H1 . The line widths are narrow (taking into account the laser bandwidth), indicating 
long lived excited states. 

These two protonated PANHs absorb in the visible or near UV spectral region, this 
absorpti on region be i ng in agreement with the absorption recorded in water solution at 
low pH (|Rvan et al.l . 119971) . This is the same absorption region as previou sly recorded 


201Oa) 

, anthracene 

flat a et al. 

2010b 

: Garkusha et al.. 

2011 ) 

. phenanthrene (Alata. 2012: 

Garkusha et al. 

2011 

), fluorine 

Alata et al.. 2012) 

, tetracene 

(Alata et al. 

2010 b). Dvrene 

jAlata 

. 2012: G 

arkusha et al.. 

201 ll 1 

4ardv et al. 

20 H 

l) and coronene ( 

Garkusha et al.. 

2011:1 

lice et al 

, 201 

4). 


Protonation of aromatic heterocycles shifts the absorption spectra to lower energy, but 
to a smaller extent than for the fully carbonated analogues: for example acridineH’*' is 
red shifted by 0.4 eV from acridine, while anthraceneH+ is red shifted by 0.9 eV from 
neutral anthracene (see Table [ij. This results in a blue shift of the transition of protonated 
acridine and phenanthridine compared to protonated anthracene and phenanthrene. The 
reyerse was observed in the case of radical cations: ionised qui noline and isoquino line have 


their transitions red shifted compared to ionised naphthalene (jDrvza et al 


2012 ). 


4. Astrophysical Implications 

In this study we have reported the electronic spectra for two protonated three-ringed PANH 
cations, derived both experimentally and theoretically. These are the first such experimen¬ 
tal spectra, and allow us to draw several astrophysically-relevant conclusions for H+PANH 


fuse interstellar band (DIB) carriers ( 

-in'--'-'- -- --- 

Salama et al.. 

1996 

Le Page et al. 

1997 

Snow et al. 

1998; 

Pathak & Sarre. 

2008 

Hammonds et al. 

200 < 

y and this study was designed to deter- 


mine if protonated PANHs could be DIBs carriers by comparing the laboratory gas phase 


absorptions at low temperature to astronomical data (IHobbs et al 


2008 


20091) . 


The first important conclusion is that H+PANH ions represent good candidates for the 
source of the DIBs. The substitution of a N atom into the H+PAH skeleton introduces 
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Table 2. Comparison of observed DIB wavelengths and equivalent widths with experimen¬ 
tal AcridineH'*' wavelengths (in nm). 


DIB 

wavelength 

442.883^ /442.819^ /442.888®"® 
437.173V436.386^ 


425.90D 


417.65®“® 


FWHM/A 

22.5^’V12.3® 

1.03/0.46^’® 


1 . 05 ^’® 


23.3® 


Experimental AcH^ 


wavelength 

intensity“ 

Assignment' 

442.24 

100 

OO 

O 

439.54 

15 

l2 

434.86 

170 

lOj 

433.38 

30 


432.74 

85 

3? 

432.18 

15 


431.42 

30 


431.04 

130 

18j 

430.30 

55 


429.52 

95 


428.9 

20 


428.3 

65 


427.62 

150 

log 

427.14 

25 


426.22 

20 


425.62 

55 


424.32 

20 


423.98 

110 

10jl8g 

423.28 

25 


423.02 

50 


422.46 

55 


420.56 

60 

log 

420.24 

55 


418.94 

25 


417.62 

50 

10gl8g 

417.16 

35 


416.24 

35 


415.46 

60 


414.84 

65 


413.83 

85 

10jl8g 

413.5 

45 


412.72 

35 


411.08 

25 



Re ferences. (1) 


Hobbs et al. 


f4 HKrelowski et al.l (119951 1: (5) 


1200911: f2l 


Hobbs et al. 


Tiiairisg et al. 


1 200811 : (3) I.Tenniskens fc DesertI (|l99J); 


Notes. The intensity of the vibronic bands relative to the origin band are given (as percent¬ 
ages), but the absolute values should be used for guidance only, as there are possible saturation 
effects in the experimental method which are not well controlled. The notation Vq corresponds 
to a transition from the vibrationless level in the ground electronic state, So, to a level with n 
quanta of the vibrational mode v in the Si electronic excited state. At the low temperature of the 
experiment, only the vibrationless level is populated in the So ground state. 
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a blue shift in the UV-visible electronic spectrum, resulting in bands at higher energies 
than those observed fo r H+PAHs. It has previo usly been shown tha t both N-substitution 
(|Hudgins et al.l . l2005h and protonation of PAHs (jAlvaro Galue et all 1201011 introduce blue 
shifts in the IR emission features, giving a better fit to the UIR features attributed to 
PAHs. A direct compariso n of our H+PANH spe ctra with H+PAH spectra of anthracene- 
9H and phenanthene-9H (|Garkusha et al.l . 120111) reveals blueshifts of ^ 11 and 133 nm, 
respectively, due to the presence of N in the ion. The destabilising effect of N on the ion in 
its excited state is particularly marked for phenanthridine. For the small PANH investigated 
here, this blue shift results in the majority of their absorption features occurring at higher 
energies than those observed for DIBs, but this does not preclude the bands of larger, more 
astrophysically-relevant PANHs falling within the DIBs frequency range. 

Importantly, the absorption lines in our spectra have FWHM of approximately 10 cm“^ 
15 - 20 A), which (upon correction for the ^ 8 cm“^ bandwidth of the laser) are 
comparable to observed DIB line widths of 2 - 3 cm“^. Our data are directly comparable 
to all of the DIBs in the 445 - 410 nm range, two of which have FWHM of ~ 23 A (> 
100 cm“^) and three of which are on the order of 1 A (2 - 5 cm“^, see Tabled]). The narrow 
bands of these H+PANH species are in di rect contrast wi t h tho se of PAH+ radical cations, 
which are much wider (~ 20 - 30 cm“^, iHudeins et al.l (120051) ') due to their intrinsically 


shorter excited state lifetimes. Due to the closed-shell nature of the H+PANH, they are also 
more photochemically stable than open-shell radical cation species, and therefore would be 
expected to be more stable in, for example, PDRs. It should be noted here that neutral 
PAHs (wh ich are themselves clo sed-shell species) also exhibit narrow absorption lines of 


3 cm ^ (Tan & Salama. 


2005^ and are likely to be major contributors to the weak DIBs. 


The first electronic transition that we observe in the spectrum of protonated acridine is 
centred at 4422.4 A and has a bandwidth-corrected FWHM of < 1 5 A. The strongest 


broad DIB is centred at 4428.2 A and has a FWHM of 22.50 A ( Hobbs et al 


2008 


■lenniskens fc DesertI [l994). The AcH'*' band is shifted by ^ 6 A compared to this un¬ 


usually broad DIB, and thus we conclude that it is not responsible for the interstellar 
feature. The other strong bands in the AcH+ spectrum do not match known DIBs either 
(see Table HI). Moreover, the poor correlation between DIBs precludes the contribution of 
molecular systems with high vibrational activity to these bands. Although the protonated 
PANHs studied here present fewer active modes than their fully carbonated analogues, they 
do show vibrational progressions. The same behaviour seems to exist i n small PANH rad i- 


cal cations, which present long vibrational progressions in the visible (iDrvza et al 


2012 ). 


Nonetheless, the progressions that we observe in our small H+PANHs would likely become 
less pronounced in larger species as they derive from vibrational modes which would be 
damped by the larger aromatic structure. Additionally, PANHs bend less than the corre¬ 
sponding PAHs, therefore there is more chance of observing isolated bands. Protonation 
of the studied PANHs occurs preferentially on the N atom, and both species contain N 
substituted in the PAH exoskeleton. The most favourable protonation orientation was thus 
H'*' in the plane of the molecule, maintaining the planar nature of the species upon ioni¬ 
sation. This would not necessarily be the case for endoskeletal N as (i) the presence of N 
already breaks the planar nature of the molecule and (ii) if protonation occurs on the N, 
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it would occur above or below the planar ring structure. The resulting non-planar struc¬ 
ture would most likely have additional vibration in electronically excited states, resulting 
in vibrational progressions in the electronic spectra. Thus, exoskeletal H+PANH probably 
represent better candidates for the DIBs. 

However, as discussed above, it is important to remember that small H+PAH and 
H+PANH molecules will not be stable in interstellar environments such as PDRs or dif¬ 
fuse HI clouds, and these experiments must be extended to species with larger skeletons. 
Previo us studies suggest tha t the 6.2 ^m band is best reproduced by PAHs with 60 - 90 C 


atoms ( Schutte et al 


199311 and with 2 - 3 N atoms (jHudgins et al.l . 1200511 . Such species 


may be photostable and their electronic transitions should be more to the red, where there 
is a larger number of DIB s. As mentioned in § [TJ the interstellar abundance of N is ^ 0.25 
that of C (ISoitzerl . 1197811 . Thus, from a purely statistical point of view, all PAHs large 
enough to be stable in the ISM (approximately > 50 C) should contain more than ^ 10 N 
atoms. Any deviation from these statistical ab undances there f ore re veals evidence of selec¬ 
tivity in the formation routes to such species. iHudgins et al.l (|2005ll conclude that “most” 
PAH species contributing to the UIR features contain nitrogen atoms, with an estimated 
> 1.2 % of interstellar N contained in PANHs . When compa red to the fraction of elemental 
C locked up in PAHs (approximately 3.5 % (jTielensl . 1200811 1. the relative abundance of N 
appears slightly more than that expected from the statistical N:C ratio, but this could be 
due to the difference in derivation met hods for PAH and PA NH abundances. 

On the basis of the results of IR ([Hudgins et al.l . 1200511 and UV-visible (this work) 
studies into PANHs, it is clear that such species are good candidates for inclusion in the 
“PAH hypothesis”. Although the spectra of the small H+PANHs studied in this work do 
not, themselves, correspond to observed DIBs, this preliminary study points to the need for 
further examination of this class of molecular ion. In particular, the effects of larger PANHs, 
substitution of larger numbers of N (or other heteroatoms), and endoskeletal heteroatom 
substitution on the UV-visible and IR spectra of protonated PANH molecules should be 
investigated. 
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Appendix A: Structures of acridine and phenanthridine 


8 9 1 



Fig. A.l. Numbering of atoms in acridine and phenanthridine used in this work (note that 
in these molecules, positions 9 and 10 are equivalent). 
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Table A.l. AcridineH+ optimised geometry for the ground (So) and first excited (Si) 
states, given in Cartesian coordinates. The columns are: the atom type; the x coordinate; 
the y coordinate; the z coordinate. All coordinates are in Angstroms. To visualise the 
geometry of the molecule, copy this table into quantum chemistry visualisation software 
(e.g. molden, molekel, gausview, gabedit, molview.) 


So 


c 

-.15633932804028E-01 

-.10862080610293E-03 

0.12980651323101E-01 

c 

0.21604311219701E-02 

-.19089158226843E-03 

0.14432770068046E-f01 

c 

0.11972021055134E+01 

0.89803276031878E-05 

0.21430993948263E-f01 

c 

0.24199216753256E+01 

0.23481628077462E-03 

0.14169482086920E-f01 

c 

0.24217097002961E+01 

0.21387517383260E-03 

-.22952726592123E-01 

c 

0.11637565441991E+01 

0.63721063504643E-04 

-.70947996359481E-f00 

c 

0.36588827993023E+01 

0.89670036868298E-04 

-.69413240849277E-f00 

c 

0.48768607333507E+01 

0.99671491006516E-04 

0.11281256027584E-01 

c 

0.48385023831464E+01 

0.11557640267375E-03 

0.14506720773022E-f01 

N 

0.36204394100487E+01 

0.26148585104960E-03 

0.20629420747075E-f01 

c 

0.61534664832050E+01 

0.39462739087322E-04 

-.63990465747766E-f00 

c 

0.73122559166118E+01 

-.85881352283063E-04 

0.11515796681357E-f00 

c 

0.60405002971524E+01 

-.34451895911368E-04 

0.22106321594476E-f01 

H 

0.12086979465670E+01 

0.27318489630465E-05 

0.32385711719424E-f01 

H 

-.94338837841307E+00 

-.43462856205420E-03 

0.19943046581554E-f01 

H 

-.97644198045926E+00 

-.12037658467841E-03 

-.51018374593142E-f00 

H 

0.11594710950112E+01 

0.92812616426702E-04 

-.18039781566228E-f01 

H 

0.36741793614736E+01 

-.14223620255338E-03 

-.17912511773748E-f01 

H 

0.59984611522114E+01 

-.24205575965616E-04 

0.33053573920256E-f01 

H 

0.82872752117088E+01 

-.10482312704436E-03 

-.38101766278469E-f00 

C 

0.72545895882594E+01 

-.16468849998998E-03 

0.15444017023691E-f01 

H 

0.61882703893466E+01 

0.12746708578654E-03 

-.17338576901536E-f01 

H 

0.81844057585583E+01 

-.36769430733588E-03 

0.21215810031686E-f01 

H 

0.36061730178587E+01 

0.47822758248642E-03 

0.30861488381914E-f01 

Si 

c 

-.70929314279567E-01 

-.31541786688268E-03 

0.44608275154867E-01 

c 

-.58897373717539E-01 

-.13395713442835E-03 

0.14283873520112E-f01 

c 

0.11957023956751E+01 

0.15274479310515E-03 

0.21215288821553E-f01 

c 

0.23990322283165E+01 

0.34800353268296E-03 

0.14062986942649E-f01 

c 

0.24044599532987E+01 

0.40680857123459E-03 

-.25754771470108E-01 

c 

0.11674457458772E+01 

0.24263938559302E-03 

-.68146830848166E-f00 

c 

0.36594857191751E+01 

0.15905476910033E-03 

-.73732905055092E-f00 

c 

0.48941828905466E+01 

-.10639780193868E-03 

0.89621262121113E-02 

c 

0.48596781329065E+01 

-. 13721149479062E-03 

0.14406105782956E-f01 

N 

0.36202263549778E+01 

0.86955555872995E-04 

0.20781472817843E-f01 

c 

0.61489971947025E+01 

-.10031790764995E-03 

-.61200844532944E-f00 

c 

0.73666490567941E+01 

-.12154090264001E-03 

0.14831140098671E-f00 

c 

0.60425997832978E+01 

-. 1065654346971 lE-03 

0.21891108230731E-f01 

H 

0.12096382514204E+01 

0.25144041827081E-03 

0.32181100284208E-f01 

H 

-.98921241395055E+00 

-.43914071065955E-04 

0.20038776168137E-f01 

H 

-.10140993189400E+01 

-.11266848263259E-02 

-.50972845019095E-f00 

H 

0.11482133821656E+01 

0.76182120393842E-03 

-.17777117058225E-f01 

H 

0.36747793130142E+01 

0.15942119712985E-03 

-.18340692796154E-f01 

H 

0.59980989482313E+01 

-.11221401574502E-03 

0.32848772163235E-f01 

H 

0.83249074127402E+01 

-.30888517421707E-03 

-.37951423721326E-f00 
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Table A.2. PhenanthridineH+ optimised geometry for the ground (Sq) and first excited 
(Si) states, given in Cartesian coordinates. The columns are: the atom type; the x coordi¬ 
nate; the y coordinate; the z coordinate. All coordinates are in Angstroms. To visualise the 
geometry of the molecule, copy this table into quantum chemistry visualisation software 
(e.g. molden, molekel, gausview, gabedit, molview.) 


So 


c 

-.29111195181980E-b01 

-.65146489862514E-b00 

-.61844760890911E-04 

c 

-.29387908722540E-b01 

-.20410167460978E-b01 

-.82569277299353E-05 

c 

-.16588107991563E-b01 

0.77733855689781E-02 

0.75458641293230E-04 

c 

-.17219813446285E-b01 

-.27709323082638E-b01 

-.37199921828215E-04 

c 

-.41886029648147E-b00 

-.69987309161869E-b00 

0.10899941288753E-03 

c 

-.49024121706251E-b00 

-.21199493808559E-b01 

-.99735744507177E-05 

c 

-.51015052010200E-b00 

0.21191676941569E-b01 

0.23906139786150E-04 

c 

0.75985855123720E-b00 

0.14869008357507E-b01 

0.45563123431349E-04 

c 

0.82224601841610E-b00 

0.51914718321358E-01 

0.86731052149833E-04 

H 

-.38388671239447E-b01 

-.68177206849797E-01 

-.29742467374202E-03 

H 

-.38969462560363E-b01 

-.25683393257787E-b01 

0.15197117974559E-03 

H 

-.17502597060224E-b01 

-.38647712115249E-b01 

-.50293301186261E-04 

H 

0.42725694029513E-b00 

-.27127154355780E-b01 

-.12097218655790E-03 

H 

-.62271962043965E-b00 

0.32076991915796E-b01 

-.14862031904349E-04 

C 

0.19418691452044E-b01 

0.22795093636996E-b01 

0.24140782153817E-04 

c 

0.21063891194358E-b01 

-.55155814726827E-b00 

0.95493200063119E-04 

c 

0.32588091784574E-b01 

0.23370986681866E-b00 

-.13841075563305E-04 

c 

0.31814283155797E-b01 

0.16504812920342E-b01 

-.66284754156874E-04 

H 

0.18625496555808E-b01 

0.33716789555816E-b01 

0.11820081928348E-03 

H 

0.22046994851041E-b01 

-.16397107772172E-b01 

0.21478512278879E-03 

H 

0.42392124513326E-b01 

-.25305100290515E-b00 

-.93056156862188E-04 

H 

0.40988144559975E-b01 

0.22465926464354E-b01 

-.22937380057018E-03 

N 

-.16236142038085E-b01 

0.13967651401319E-b01 

0.41276847684163E-04 

H 

-.25207748385068E-b01 

0.18893664425044E-b01 

0.19856844428137E-04 

Si 

c 

-.28919722947996E-b01 

-.65278718807287E-b00 

0.33414688738087E-04 

c 

-.29004845031867E-b01 

-.20743437122434E-b01 

-.25288045960361E-04 

c 

-.16856282174091E-b01 

0.51021196541972E-01 

0.93990452289923E-04 

c 

-.16975565658929E-b01 

-.27971634118656E-b01 

-.51094140371324E-04 

c 

-.43113733728882E-b00 

-.66970374682184E-b00 

0.94905579705874E-04 

c 

-.48109651644802E-b00 

-.20981763650497E-b01 

0.10747059410705E-04 

c 

-.53502759597706E-b00 

0.21924497474230E-b01 

-.20002232469464E-04 

c 

0.73493449624166E-b00 

0.15208906115057E-b01 

-.14872482318922E-04 

c 

0.80956248952303E-b00 

0.80419937319906E-01 

0.43867431595168E-04 

H 

-.38417639939713E-b01 

-.10589960320023E-b00 

0.24604203282430E-04 

H 

-.38625961765093E-b01 

-.25968695105859E-b01 

-.21430470646605E-04 

H 

-.17040080213666E-b01 

-.38902513999984E-b01 

-.11837063774503E-03 

H 

0.45192504967185E-b00 

-.26674086478333E-b01 

-.21868947529171E-04 

H 

-.63965656567796E-b00 

0.32799929839071E-b01 

-.66438530599537E-04 

C 

0.19303418180185E-b01 

0.22476705272306E-b01 

-.50108855842193E-04 

c 

0.20903107675528E-b01 

-.54987302878548E-b00 

0.37158316542311E-04 

14:1 

0.32945230481049E-b01 

0.19861737786567E-b00 

0.20321455353585E-04 

c 

0.32072629909813E-b01 

0.15900369029004E-b01 

-.10390058308755E-04 

H 

0.18968154620148E-b01 

0.33434241873007E-b01 

-.11722732237866E-03 

H 

0.21566285807046E-b01 

-.16413683816682E-b01 

0.23874601091976E-04 
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Appendix B: Mass spectra 






Electron impact ionisation of acridine 

-f - 

140 



T 


160 


Fig. B.l. AcridineH+ photofragmentation mass spectrum and comparison with collision 
induced fragmentation and electron impact spectrum of neutral acridine. The first panel 
shows the experimental mass spectra: the mass spectrum without fragmentation laser is 
shown in red (and is multiplied by -1 to facilitate comparison), and the mass spectrum 
recorded with the fragmentation laser on is shown in black. The black spectrum includes 
the newly formed fragments around m/z 152 (major fragment) and around m/z 127 (minor 
fragments). The violet curve presents a zoom of the region around m/z 152, presented as a 
difference spectrum of the recorded spectra with laser on - laser off. This photofragmenta¬ 
tion spectrum is different from the mass spectrum obtained by collision induced dissociation 
of AcH+ (plotted in black), in which the fragments m/z 150 and 151 have much higher 
intensity. It also differs from the mass spectrum obtained by electron impact ionisation of 
neutral acridine (plotted in red) which peaks at m/z 151. 
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Mass 


Fig. B.2. PhenanthridineH’*' photofragmentation mass spectrum: the red trace corresponds 
to the mass spectrum recorded without the fragmentation laser, the black trace is the mass 
spectrum recorded with the fragmentation laser on and shows the appearance of fragments 
around m/z 152 (major fragments) and around m/z 127 (minor fragments). One spectrum 
is positive and the other negative just to present clearly the difference between the two 
spectra. 
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Table C.l. PhenanthridineH+ vibronic bands: only the bands with large intensities are 
tabulated. Above 26500 cm“^, the assignments are only tentative. 


Wavelength (nm) 

Wavenumbers (cm ^) 

Energy above origin (cm 

Intensity 

Assignment 

397.28 

25171 

0 

m 

o§ 

393.48 

25414 

243 

m 


389.8 

25654 

483 

m 

5§ & 15j 

389.04 

25704 

533 

str 

17j 

386.16 

25896 

725 

m 

5g & 15j5j 

385.44 

25944 

773 

str 

5jl7j 

382.68 

26131 

960 

m 

5^ 

381.84 

26189 

1018 

str 

5§17j 

381.16 

26236 

1065 

str 

17o 

379.64 

26341 

1170 

m 


378.36 

26430 

1259 

m 

5§17j 

377.68 

26477 

1306 

str 

5jl7i 

376.72 

26545 

1374 

str 


376.2 

26582 

1411 

str 


374.92 

26672 

1501 

m 


374.28 

26718 

1547 

str 

5§17g 

373.56 

26769 

1598 

str 

17g 

372.12 

26873 

1702 

str 


371.68 

26905 

1734 

m 


370.88 

26963 

1792 

m 


369.88 

27036 

1865 

str 


369.4 

27071 

1900 

str 


368.84 

27112 

1941 

str 


367.64 

27201 

2030 

m 


366.2 

27307 

2136 

str 

17^ 

362.24 

27606 

2435 

str 



Appendix C: PhenanthridineH^ vibronic bands 
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